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Abstract: 

The decay B — > K*{-^ Ktt)U offers great opportunities to explore the physics at and above 
the electroweak scale by means of an angular analysis. We investigate the physics potential 
of the seven CP asymmetries plus the asymmetry in the rate, working at low dilepton mass 
using QCD factorization at next-to leading order (NLO). The b — * s CP asymmetries are 
doubly Cabibbo-suppressed < 1% in the Standard Model and its extensions where the 
CKM matrix is the only source of CP violation. Three CP asymmetries are T-odd, and 
can be 0(1) in the presence of non-standard CP violation. The T-even asymmetries can 
reach 0(0.1), limited by the small strong phases in the large recoil region. We furthermore 
point out an easy way to measure CP phases from time-integrated, untagged B^jB^ — ► 
K*(-^ K°ir°)ll and B S ,B S — > (p(— > K + K~)ll decays. Analyses of these CP asymmetries 
can rule out, or further support the minimal description of CP violation through the CKM 
mechanism. Experimental studies are promising for (super) flavor factories and at hadron 
colliders. 
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1. Introduction 



The quest for physics beyond the Standard Model (SM) is of highest priority at current 
and future flavor facilities [1]. A promising direction is to look for CP symmetry breaking 
effects that cannot be accounted for with the SM's very predictive CKM-mechanism of CP 
and flavor violation for the quarks. Rare b — > s transitions are such sensitive probes, since 
all CP violating effects in the SM receive a universal suppression of order 10~ 2 from the 
CKM matrix elements. Ongoing experimental efforts include those with hadronic decays 
such as B S ,B S — > J/^cj) at the Tevatron [2] and b — > sss induced penguin modes, see, e.g., 
[3] for recent data. 

We focus here on exclusive semileptonic decays induced by b — > sll, I = e, fj,, to explore 
the borders of the SM. The decays B — ► Kll and B — > K*ll have already been observed, 
with their rate being in agreement with the SM. Investigations of more involved observables 
such as dilepton mass spectra, lepton angle distributions and dimuon to dielectron ratios 
are currently being pursued [4, 5, 6]. The decay B — > K*ll with subsequent K* — > Kir 
allows to extract further physics information through an angular analysis of the Kirll final 
state [7], especially on the chirality content [8, 9, 10]. 

Here we study the eight CP asymmetries in B — ► K*{-^ Kir)ll decays, the one in 
the decay rate plus seven more requiring angular information. We calculate the matrix 
elements to next-to-leading order (NLO) in the strong coupling and to lowest order (LO) 
in the 1/E expansion [11], where E denotes the energy of the emitted K* in B rest frame. 
We work within QCD factorization (QCDF), which has been applied to B — ► K*ll decays 
[12, 13], for analyses in soft-collinear effective theory (SCET), see [14]. Previous works on 
CP asymmetries in the angular distributions [7, 15, 16] employed naive factorization, see 
[17] for a PQCD study. 

By means of the large energy limit, the QCDF framework used in our analysis holds 
for small dilepton invariant masses, and predicts small strong phases: Lowest order quark 
loops are either close to or below threshold, hence real (charm quarks), or CKM suppressed 
or induced by small penguin contributions (the lighter quarks). Other sources of strong 
phases are from subleading spectator effects or higher order in a s . 

This highlights T-odd CP asymmetries, three of which are accessible with the angular 
analysis. (The T-transformation reverses the sign of all particle momenta and spins.) The 
important feature is that the T-odd asymmetries are oc cosAgsinAvy, where As and 
Aw denote the differences of strong and weak phases, respectively. While the T-even 
CP asymmetries oc sin As sin A vanish for small strong phases, the T-odd asymmetries 
exhibit maximal sensitivity to CP violation in this limit. 

Time-dependent CP asymmetries in Bd,Bd — ► K*(— > K°tt )11 have also been con- 
sidered in [16]. These measurements require large amounts of data. We suggest here to 
use the CP-odd property of four of the asymmetries to extract them from an untagged 
and time-integrated data set. It has been known that such data sets are useful to access 
CP violation in angular distributions, e.g., [18, 19]. We work out the corresponding CP- 
sensitive observables in B^B^ — ► K*{— ► K°tt°)11 and B S ,B S — > <p{— > K + K~)ll decays. 
The K + K~ll final state is the same as the one of the B S ,B S — >■ </>(—>■ K + K~)J/^>{-^ 11) 
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decays, which are already under experimental study including angular analysis [2]. 

In Section 2 we review the B — > K*(—> Kir)ll angular distributions. CP asymmetries 
and possibilities of their measurement from double- and single-differential distributions are 
given in Section 3. Prospects for B d ,B d -► K*(^ K°ir°)ll and B S ,B S -► </>(-> K + R-)ll 
decays without tagging are examined in Section 4. SM predictions and theoretical un- 
certainties are presented in Section 5. In Section 6 we model-independently investigate 
the impact of New Physics (NP) on the CP asymmetries. To do so, we work out con- 
straints from other rare decay data, B — > X s ll,X s j, the forward-backward asymmetry in 
B — > K*ll and time-dependent CP asymmetries in B d ,B d — > ► K°tt )^. We sum- 

marize in Section 7. Various appendices A-F contain details of the calculation of the CP 
asymmetries. 

2. Full Angular Distribution 

In this section we review the angular distribution of the exclusive decay B — ► K*{-^ Kir)ll 
and its conjugate decay. Throughout this work we use B = (bq) for q = u,d, B s = (bs) 
and K,K* = (sq). We illustrate the kinematics for neutral mesons decaying to charged 
particles. Charged -B-decays can be treated analogously. We discuss meson mixing effects 
and B s — ► <j>{— > K + K~)ll decays in Section 4. 

The full angular distribution of the decay B° -► K--k + )11 can be written in the 

limit of an on-shell decaying K* resonance as a 4-differential distribution [7, 9] 

f/Jj ^J{q 2 ,9 h e K *^), (2.1) 



dq 2 d cos 9id cos 6k* d(f) 8tt" 

where the lepton spins have been summed over. Here, q 2 is the dilepton invariant mass 
squared, that is, g M is the sum of and pf_, the four momenta of the positively and 
negatively charged lepton, respectively. Furthermore, 9i is defined as the angle between 
the negatively charged lepton and the B in the dilepton center of mass system (c.m.s.) and 
9k* is the angle between the Kaon and the B in the (K~tt + ) c.m.s.. We denote by Pi the 
three momentum vector of particle i in the B rest frame. Then, (j) is given by the angle 
between p K ~ x p w + and p ; x p;+, i.e., the angle between the normals of the (K~tt + ) and 
planes. The full kinematically accessible phase space is bounded by 

Am 2 < q 2 < {M B - M K *) 2 , -1<cos0j<1, -1 < cos 9 K * < 1, < <p < 2vr, (2.2) 

where mi, Mb and M^-* denote the mass of the lepton, B meson and the K*, respectively. 
For an off-resonance B — > K-zr// study, see [20] . 

The dependence of the decay distribution (2.1) on the angles 9i, 9k* and 4> can be 
made explicit as 



J(q 2 , 9i, 9 K * , <t>) = J{ sin 2 9 K * + J{ cos 2 0^* + ( Jf sin 2 + J| cos 2 6>x» ) cos 29, 



+ J3 sin 2 sin 2 9\ cos 2(f) + J4 sin 20R-* sin 20; cos <fi + J5 sin 20^-* sin 6*/ cos 
+ Jq sin 2 0^* cos 9\ + J7 sin 20^* sin 9\ sin 

+ J$ sin 20R-* sin 26*/ sin <p + Jg sin 2 0^* sin 2 9\ sin 20, (2-3) 
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where the coefficients = J- a ^ (q 2 ) for i = 1, . . . , 9 and a = s,c are functions of the dilep- 
ton mass. In the following we suppress the g 2 -dependence for brevity also in expressions 
derived from the . 

The angular coefficients can be expressed through the K* transversity amplitudes 
Ai(q 2 ) with i = {_L, ||,0,t}, see Appendix A. Note that not all the are independent, 
for example, for vanishing lepton masses 

3Jf = J|, J{ = -J 2 C . (2.4) 

Furthermore, in the absence of right-handed currents J3 = Jg = up to power corrections. 

The corresponding distribution of the CP conjugated decay B° — > K*°(^ K + ir~)U 
can be written as 

d 4 T 3 -, 2 



—j{ q \e h e K ,A)- (2.5) 



dq 2 d cos 6 id cos 9 k* d(p 8tt 

Here, 9k* denotes the angle between the Kaon and the B meson in the (K + tt~) c.m.s.. 
The definiton of 9\ is identical for both B and B decays. The angle <j) for B decays is given 
by the angle between p K + x p^ and p ; x p ; + . Therefore, in the limit of unbroken CP, the 
distributions for B and B mesons are equal under the combined transformations 9i — > 6*/ — 7r 
and <fi — > —(j). The function J is hence obtained from J in (2.3) by the replacements 

J h2,3,4J Ji%3,4J> ^5,6,8,9 ~> "^5,6,8,9, (2-6) 

where equals with all weak phases being conjugated [7]. 

With its rich multi-dimensional structure, the angular distributions in (2.1) and (2.5) 
have sensitivity to various effects modifying the SM, such as CP violation beyond CKM 
and/or right-handed currents. Given sufficient data, all and can in principle be 
completely measured from the full angular distribution in all three angles 9[, 9k* and (f>. 

The familiar dilepton invariant mass spectrum for B — ► K*ll decays can be recovered 
after integration over all angles as 

^2=«/i-y, where Ji, 2 = 2Jf >2 + J^. (2.7) 

The (normalized) forward-backward asymmetry ^4fb is given after full 4> and 9k* integra- 
tion as 1 



A FB (q 2 



d cos 9\ 



d 2 T 



dr 



dq 2 d cos 9i / dq 2 



dq 2 ' 



(2.8) 



By dT /dq 2 and Afb(q 2 ) we refer to the corresponding spectra of the CP conjugated decays. 



1 Since we define the lepton angle 8i with respect to the I , our definiton of the forward-backward 
asymmetry (2.8) differs from the one in other works using the l + , e.g., [5, 6, 12, 13, 21], by a global sign. 
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3. CP Asymmetries 

CP violating effects in the angular distribution are signaled by non-vanishing differences 
between the (g 2 -dependent) angular coefficients 

Aj\ a) =Aj\ a \q 2 ) = 4 a) -j\ a) . (3.1) 

Of particular importance are the asymmetries related to the coefficients ^7,8,9- These 
are odd under <fi — > —<f>, and hence induce T-odd asymmetries AJ7 8 9 which are not sup- 
pressed by small strong phases as predicted from QCDF. 

The CP asymmetry in the dilepton mass distribution is commonly defined as, see (2.7), 



a (n 2,_ d(T-T) Mr + r) 1 



AJi - ^ 



, iVr = N r (q 2 ) = ^fi. (3.2) 



Following [7], we define in addition to Aqp seven normalized CP asymmetries as 

oA j 9A j. 

Mq2) = lvf for * = 3 ' 6 ' 9 ' A ^ q2) = -lvf for * = 4 > 5 >7,8. ( 3 - 3 ) 

Note that up to differences in the normalization A 6 equals the forward-backward CP asym- 
metry Apg advocated to search for non-standard CP violation in the decay B — > K*ll 
[22, 23], see (2.8). For ^-integrated quantities we introduce the notation 

(X)= r™dq 2 X(q% (3.4) 

J a 2 . 

where the integration from q^ in to q^ax should be in the low dilepton mass region in order 
to use the 1/E expansion of QCD for theory predictions, see Appendix C. We then define 
the normalized g 2 -integrated CP asymmetries as 

(A i ) = 2 { -^1 fori = 3, 6, 9, (A? > = for i = 4, 5, 7, 8, (3.5) 

where the numerator and the denominator are integrated with the same q 2 cuts. 

The CP asymmetries (Ai) (i = 3, 6, 9) can, for example, be extracted from the double- 
differential distribution in d\ and (f>, 

, ^ , , = ^-{ (Ji) + (J2) cos 20, + 2 ( J 3 > sin 2 0, cos 20 
a cos 0; a<p 47r 

+ 2 ( J 6 ) cos 0, + 2 ( J 9 ) sin 2 0j sin 20} , (3.6) 
which is obtained from integrating (2.1) over Ok*- After full ^^-integration follows 

" M ' : -^((Ji)-^# + ^(J3)cos20 + l(J 9 )sin20l, (3.7) 



# 27r ^ 3 3 "3 

showing that (AJ9) can be found from d^T + f ) /dcj), whereas (AJ3) can be obtained from 
<i(r — f ) /d(f), with (AJi) — (AJ2) /3 from Acp without angular study, see (3.2). 
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The measurement of the CP asymmetries {Af) (i = 4,5,7,8) requires binning into 

cos Ok* as 



p^A \f.n dma e K . *m tM ,3.8, 

acosOidcp IJ J-ij d cos 9k* d cos t>i d<p 

= — { (J4) sin26*i cos <fi + (^5) sin #2 cos <fi + (J7) sin^ sin^> + ( Jg) sin 26*/ sin</>}. 

From here follows upon full ^-integration 

d(A e \ 1 



dcp 4 



{(J 5 )cos0 + (J 7 )sin0}. (3.9) 



We learn that (AJ5) can be extracted from d(Ag Kt + Aq k ^ /dcf) whereas (AJ7) can be 
obtained from d(Ao Kt — Aq k , ) /g?0. The double asymmetry in 9k* and 6>/, 



-ir-/: 



d 2 (A ) 2 
dcos °i a ink = T~ cos0+ (J 8 ) sin</>} , (3.10) 
d cos 9\ dep 37r 



allows to obtain (AJ 4 ) from d (^4e K , ,e t ~ ^-0 K *,6i) /d<t>, whereas (AJ 8 ) can be extracted 
from d (Ao K , fil + A~e K » ,e t } /#■ 

The latter considerations demonstrate how the CP violating angular coefficients A J; 
for i = 3, 4, 5, 7, 8, 9 can be extracted from distributions in the angle 4>. The quantity A Jg 
can be measured easiest from the cos ^-distribution, i.e., by adding the (numerators of 
the) forward-backward asymmetries Afb and ^4fb, see (2.8). Note that only A3, A§ and 
Ag can be obtained from a genuinely single-differential distribution. Ag is the only T-odd 
asymmetry with this property. 

Analytical expressions for all CP asymmetries at NLO in terms of the short distance 
coefficients from the electroweak Hamiltonian in Appendix B are presented in Appendix D. 
In particular, we include NLO a s -corrections thus present in this work the first analyses of 
CP asymmetries in the B — ► K*{-^ Kir)ll angular distributions at this order. 

Corrections from Bj — B^ mixing to the time integrated CP asymmetries in flavor 
specific (self-tagging) final states are of the order \ A$ L \ < O(10 -3 ) [24] and can be neglected. 
Here, A§ L denotes the semileptonic asymmetry into wrong sign leptons in the S^-system. 

4. CP Asymmetries without Tagging 

With the Jj, i = 5, 6, 8, 9 being odd under CP, see (2.6), the corresponding CP asymmetries 
can be extracted from dT+dT, i.e., without identifying the flavor of the initial b quark. This 
feature is very useful for B d ,B d — »■ K*(-^ K°ir°)ll and B S ,B S — »■ </>(—► K + K~)ll decays, 
which unlike Bd,Bd —> K*{-^ K^ir^ll or charged i?-decays, are not self-tagging. We 
focus here on i? s -decays to CP eigenstates, but the formalism equally applies to S^-decays 
after the corresponding replacements. 

Both B s and B s angular distributions are described by the angles 9i, 9k* and <j>. 
These are defined in complete analogy with B — > K*(-^ Ktt)U decays, see Section 2: 9[ 
is the angle between the negatively charged lepton and the B s /B s in the dilepton c.m.s., 
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Ok* denotes the angle between the K~ and the B s /B s in the (K~K + ) c.m.s. and (ft is the 
angle between Px x Pk+ an d Pi- x Pi+- 

To account for mixing, time-dependent transversity amplitudes need to be introduced: 

A a (t) = A(B s (t)^4>(^K + K-)Jl), A a (t) = A(B s (t)^4>(^K + K-)Jl), (4.1) 

where A a (t), (A a (t)) denotes the amplitude for a meson born at time t = as a B s , (B s ) 
decaying through the transversity amplitude a =_L,||,0 at later times t. Here we use 
for brevity A a (t) for both A^(t) and A^(t). The formulae for the unmixed transversity 
amplitudes, i.e., the ones at t = can be taken from Appendix C with the requisite 
replacements in masses and hadronic parameters and differences in the spectator effects 
given in [13] to account for the B s — ► <\> transitions. 

Untagged rates dF + dT can then be written as (a, b =_L, ||, 0) [18] 

At(t)A a (t) + Al(t)A a {t) = ^A(B S - 0(- K + K-) b ll)*A{B s - 0(- K + K-)Jl) 



x 



[(1 + %%£*£ a ) (e" rL ' + e- r «') + + r/ a £ a ) (e" rL * - e"^)] • (4.2) 



Here, 77^0 = +1 and r/j_ = — 1 are the CP eigenvalues of the final state and ^l(h) denotes 
the width of the lighter (heavier) mass eigenstate. We neglect CP violation in mixing, 
which is bounded by the semileptonic asymmetry in the i? s -system |-A| L | ~ 0(1CT 2 ) [2]. 
Furthermore, 

c A ^s -> K+K~)Jl) 

?a A(B S - <f>(-^ K+K-)Jl)(Sw - -SwY { • ) 

where {5w — * — &w) implies the conjugation of all weak phases in the denominator and &m 
denotes the phase of the B s — B s mixing. It is very small in the SM, & = 2 &rg(V t * s V tb ) . 

The CP asymmetries AJj(i) = Jj(t) — Ji(t), i = 5,6,8,9 are then obtained by taking 
the real and imaginary parts of (4.2), adding or subtracting A%*(t)A%(t) + A%*(t)A%(t) for 
k = L and k = R, and taking into account normalization factors depending on the angular 
coefficient Jj, see Appendix A. After time-integration follows from (4.2) 

POO 

/ dt(A* b (t)A a (t) + A* b (t)A a (t)) = A(B S - 0(- K + K-) b ll)*A(B s - K + K')Jl) 
Jo 

1 Ar 

X Y(l-yl) ^ + ^bttta) ~ V {VbCb + Vata)] , V = jp, (4-4) 

where T = (Tl + Tjj)/2 and the width difference Ar = — r#. 

This expression (4.4) becomes transparent if one neglects strong phases, where £ a = 
e -i(* M -2*„) j $a = arg (^(i5 s K + K-)Jl). For a ^b,rj a = -rj b we obtain 

/ dtlm(^(t)i a (t) + ^(t)4,(t)) = ^ ~\A(B a -+ K + K-) b ll)\ 

Jo L \ l - y ) 



x\A(B s - K + K-)Jl)\ • M*a - $ 6 ) - OT a sin($ M - $a - *&)] , (4.5) 

/•oo o 

/ tit E(ia^)i 2 + i^wi 2 ) = m _ 2) E - ^-yoi 2 

x [l-2/7/ cos($ M -2$ a )], (4.6) 
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where (4.5) gives the asymmetries related to J$ t g and (4.6) is needed for normalization, 
that is, gives (twice) the CP averaged decay rate. It also exhibits sensitivity to CP phases. 
The T-even asymmetries associated with J^^ vanish with no strong phases present. 

We define (g 2 -dependent) CP-odd CP asymmetries J^^ m%x as ( a = j_ j || ; o, k = L,R) 



A mix (a 2 ) = 2 J dtAJjjt) _ 

which match the CP asymmetries of the flavor-specific, unmixed decays (3.3) for y — ► 0. 

For S^-mesons, y is below 10~ 2 [24], and the untagged and time-integrated K°ir° final 
states yield the same information on CP violation as the ones with K^ir^ discussed in 
Section 3, or charged -B-decays. (For early works with y = 0, see [19]). For l? s -mesons the 
width difference is larger, y ~ O(0. 1) [24], and interference effects become observable with 
the mixing phase $m- The latter is currently under intense experimental study and only 
poorly determined to date, see, e.g., [2]. 

We refrain in this work from presenting a dedicated numerical analysis for the B S ,B S — > 
<fi(—> K + K~)ll observables A^ mix : The presumably dominant part independent of the 
width difference can be inferred from B — > K*{-^ Kir)ll decays by 51/(3). The biggest 
corrections such as those from the form factors and phase space are expected to cancel in 
the asymmetries. On the other hand, discrepancies in the CP asymmetries between Bj- 
and ^-processes at 0(y) can be attributed to the mixing parameters y and $m- 



5. Standard Model Predictions 



CP asymmetries in the decays of hadrons are in the SM solely induced by the CKM matrix. 
For the b — > s transitions under consideration here, the requisite weak phase difference 
stems from \ u = V ub V* s /V tb V t * s . Therefore, all CP asymmetries in B — > K*(-^ Kir)ll decays 
discussed here receive an overall suppression by Im[A M ] ~ f/A 2 of order 10 -2 , where A and 
f\ denote parameters of the Wolfenstein parametrization of the CKM matrix. 

We work out the SM CP asymmetries in B — ► K*{-^ Kn)ll decays in the low-g 2 region 
using QCDF at NLO in a s and leading order 1/E. Analytical expressions for the asymme- 
tries are given in Appendix D. The CP asymmetries in the SM can be obtained by setting 
the NP Wilson coefficients Cj g 10 = 0, see Appendix B for the effective Hamiltonian used. 
Details on the QCDF framework and the transversity amplitudes are given in Appendix C. 
We take the B — > K* form factors from light cone QCD sum rules (LCSR) calculations 
[25], see Appendix E. Our numerical input is compiled in Table 1. We checked that our 
findings for the branching ratio and the forward-backward asymmetry of B — ► K*ll decays 
agree for the given input with [12, 13]. Our predictions always refer to neutral £?-decays 
unless otherwise stated. 

The three main uncertainties in the asymmetries come from the form factors £y and 
the variation of the renormalization scale fib and the CKM parameters. We vary the 
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A = 0.2258lg;{j{j$ (95% C.L.) [26] 

\V cb \ = 0.0417 ±0.0013 (95% C.L.) [26] 

p= [0.108, 0.243] (95% C.L.) [26] 

fj = [0.288, 0.375] (95% C.L.) [26] 

a s (m z ) = 0.1176 ± 0.0020 [24] 
cte(mb) = 1/133 

m w = 80.403 GeV [24] 



m 



pole 



m'c 



pole 



(170.9 ± 1.8) GeV 
(4.6 ±0.1) GeV 
(1.4 ±0.2) GeV 



[27] 
[12] 



B(B -► XM) = (10.57 ± 0.15)% [24] 
t b o = (1.530 ± 0.009) ps [24] 
Tb. = (1.425 ±0.041) ps [24] 
/ B „, d = (200±30) MeV 

f Bs = (240 ± 30) MeV [28] 
Ab,+ (1.5 GeV) = (0.458 ± 0.115) GeV [13, 29] 
/f *(1 GeV) = (185 ± 10) MeV [30] 
ff* = (217 ± 5) MeV [24] 

[25] 
[25] 



"1,K 



,(1 GeV) = 0.1 ±0.07 
^.(1 GeV) = 0.1 ±0.1 



Table 1: The numerical input used in our analysis. We denote by m\, the PS mass at the factoriza- 
tion scale [if — 2 GeV. We neglect the strange quark mass throughout this work unless otherwise 
stated. The numerical input for the form factors is given in Appendix E. 



scale between and 2rrib and allow for an uncertainty of 11% and 14% for £j_ and £m, 

respectively. The CKM input is given in Table 1. For the total uncertainty estimate, all 
three sources of uncertainty are added in quadrature. 

In Figure 1 we show the T-odd CP asymmetries A® 8 and the T-even ones ^cp , A% 5 
and Aq as a function of q 2 . The various bands indicate the uncertainties due to the form 
factors, the CKM parameters, fi b an d the total uncertainty. The asymmetries A3 and Ag 
are not shown, since they vanish in the SM at lowest order in 1/E. (A small finite value is 
induced by the strange quark mass.) Hence, their leading contributions may arise as 



^ 3 ,9~Im[A n ]0(A/£)~0(lO- 3 ). 



(5.1) 



The LO predictions for the CP asymmetries are also included in Figure 1. The higher 
order a s -corrections increase the size of the CP asymmetries. For A® and A$ this happens 
because their respective LO values are suppressed by cancellations. Specifically, in the SM 



A? ~ Im[A u ]Re 



r, 



(«) 



,2 X 



(«) 



+ 



Ml 



(5.2) 



which vanishes at LO in QCDF, see (C.4), and also [15]. (Our value of A® at LO is tiny but 
finite since in the numerical analysis we do not neglect kinematical factors M^„/M B .) The 
asymmetry A$ is subject to similar cancellations, although here an additional LO term 
exists, which is, however, numerically subleading. The values of A® and Ag are therefore 
determined by the NLO a s -corrections resulting in a large fj, b uncertainty. The impact of 
the higher order terms on the T-even asymmetries is sizeable, but less pronounced. We 
discuss further details of the SM CP asymmetries in the context of the integrated CP 
asymmetries. 

We find that all (^-integrated CP asymmetries (A^) are less than O (l0~ 2 ) in the 
SM. This can be seen in Table 2, where we give the results for the two cuts {q^ n , (/max) = 
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total(1.0) ■ 










toUl(NI-O) 




















0.004 
0.002 


-0.002 
-0.004 
-0.006 
-0.008 
-0.01 











total(NLO) 










V CKM 
total(LO) 1 




Figure 1: The CP asymmetries Acpfi(q 2 ) and Af 5 ? 8 (q 2 ) in the SM in the low-g 2 region at LO 
and NLO in QCDF. The various bands show the uncertainty due to the form factors, the CKM 
parameters and fib separately, whereas the overall band indicates the total uncertainty. 



(1,6) GeV 2 (upper entries) and (1,7) GeV 2 (lower entries), respectively. The uncertain- 
ties from the form factors and jib are also shown separately. Due to the overall CKM 
factor Im[A u ], all CP asymmetries suffer from a universal 15% uncertainty related to 
that. The corresponding values for the NLO CP averaged decay rates are (iVr)/2 = 
(1.1 ±0.3) • 10~ 19 GeV and (N T )/2 = (1.3±°;f) • 1(T 19 GeV for (q 2 min , q 2 max ) = (1,6) GeV 2 
and (1, 7) GeV 2 , respectively, assuming B(K* — > Kir) = 100%. The uncertainty in the rate 
is about 25 % from the form factors, 7 % from Vt s and order one percent from 

The form factor induced uncertainty in the asymmetries depends on the amount of 
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r i+2.4 +11 +42 
J - 1 -1.6 -8 -26 
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^• D -1.4 -6 -25 


< io- 2 


7 1+2.6 
' ■ A — 1.9 

a rr+2.3 


(Af) 


o r+1.4 +7.4 +37 
°-°-2.0 -10 -53 
o i +1.3 +6 +37 
°- i -1.7 -10 -53 


n 9+0.04 
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o 1+1.4 
J -^-2.0 
o 1+1.3 
°- i -1.8 


(A 3 ,^ 


0(1) 


0(1) 


0(1) 



Table 2: SM predictions for the integrated CP asymmetries in units of 10~ 3 with the integration 
boundaries (<7 min , '/max) = (L 6), (1, 7) GeV 2 (from top to bottom). We take into account uncertain- 
ties from the form factors, the scale dependence fib and the CKM parameters, all of them added in 
quadrature for the total uncertainty. The relative uncertainties due to £m and fib are also shown 
separately. The asymmetries at LO in a s and the NLO ones for charged i?-decays are given as well, 
see text for details. ^The leading contributions (^3,9) in the SM are power counting estimates only. 

cancellations between the numerator and the decay rate in the denominator. We recall 
that we vary the two form factors within their uncertainties independently. Taking into 
account correlations would reduce the errors in the ratios, but requires control over the 
parameters of the LCSR calculation [25], which is beyond the scope of this work. Since the 
decay rate is dominated by the longitudinal K* polarization driven by £11 , see the discussion 
following (C.7), the strongest form factor uncertainty is seen in Aq being proportional to 
The other asymmetries Af, i = 4, 5, 7, 8 with numerator oc receive more efficient 
cancellations. 

As can be seen from Table 2, (-Acp)j (A^) and (A$) exhibit a massive fib dependence of 
order 50 %. The CP asymmetries A^ with i = 4, 5, 6 are not subject to the cancellations 
mentioned after (5.2) and have a smaller residual fib uncertainty below ten percent. The 
fib dependence of (Aq) of a few percent is accidentally small due to significant cancellations 
between different g 2 -regions, see the crossing of the fib bands in Aq near q 2 ~ (3 — 4) GeV 2 
in Figure 1. 

Furthermore, we study the impact of higher order contributions in QCDF on the 
{A\ D) ). The shift from LO in a s to NLO is substantial. Switching off the spectator 
interactions reduces the size of (A® 8 ) by about 10 %, and by less for the other asymmetries. 
We also give in Table 2 the NLO SM predictions for charged i?-decays. The splitting 
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observable 


SM 


data 


B(B - X sl f 

ab 


(3.15 ±0.23) • 10" 4 [34] 

(-2.818:1) • io- 2 


(3.52 ± 0.25) • 10" 4 [3] 
-0.19 ±0.23 [3,31] 


B(B ^ X s ll)\ [lfi] 
B(B — >■ X S Z/) [ >0 .04] 


(1.59 ±0.11) • 10~ 6 [35] 
(4.15 ±0.70) • 10~ 6 [21] 
< 


(1.60 ± 0.51) • 10~ 6 [33] 
(4.5 ± 1.0) • 10 -6 [24] 
-(0.76lg;gj ± 0.07) [6], also [4, 5] 


6(5 S - w) 


~ 3 • 10~ 9 


< 4.7 • 10~ 8 at 90% C.L. [32] 



Table 3: Relevant b — > S7 and b — > sll observables. "With photon energy cut E~, > 1.6 GeV. b SM 
value obtained with m s = 0.12 GeV. c Note the different lepton angle convention between [5, 6] and 
this work. 



between the CP asymmetries in neutral versus charged -B-decays is dominated by weak 
annihilation contributions from current-current operators and varies a lot in size: (Aj? 7 ) 
((Aqp)) increase by 0(30%) (a factor of two) from neutral to charged -B-decays, whereas 
(A^) decreases by ~ 1/2. The splitting for (Agg) is at the few percent level. 

The SM predictions for the untagged, time-integrated CP asymmetries (A^ D ^ mix ), i = 

5, 6, 8, 9, in B s , B s — > K + K~)ll decays (4.7) can be inferred from the (^4^) in neutral 
S^-decays, which are given in Table 2. Corrections arise from SU(3) flavor breaking, which 
is expected to be small in the ratios, from the S s -width difference at the level of ten percent 
and from spectator interactions. All these effects are subdominant with respect to the 
theoretical uncertainties of the SM predictions. 

6. Beyond the Standard Model 

This section contains the model-independent analysis of the CP asymmetries. We consider 
NP contributions to the operators ©7,9,10 which are part of the effective Hamiltonian (B.l) 
of the SM, as well as NP contributions to the chirality flipped ones O' 7 910 . We allow the 
respective NP coefficients Cf p and C- NP = C[ for i = 7, 9, 10 to vary in magnitude and 
phase, denoted by fa, within the constraints from the FCNC B-decay data. The radiative 
and semileptonic b — > s transitions are the most important ones for our analysis. The 
relevant data and SM predictions are given in Table 3. 

In our analysis the NP Wilson coefficients are leading order coefficients. All Wilson 
coefficients are understood as evaluated at the low, ^-scale. We start with a discussion of 
the experimental constraints. 

6.1 Experimental Constraints 

The radiative decays induced by b — > S7 probe the electromagnetic dipole coefficients 
C 7 and C'-j. In the NP scenarios considered here, the flipped dipole coefficient has no 
interference terms in the radiative decay rates. Hence, these observables constrain only the 
magnitude of C' 7 and not its phase. 
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We take into account the B — > X s 7 branching ratio for which we adopt the NNLO SM 
results from [34] . To account for the missing higher order calculation of the beyond-the-SM 
amplitude, we take for the theoretical uncertainty of the NP contribution twice the SM 
uncertainty. We apply the experimental constraints at 90 % C.L. We checked that the 
direct CP asymmetry in B — > X s j, e.g., [36], does not give constraints beyond those from 
the B — > X s 7 branching ratio. 

The time-dependent CP asymmetry Sk*^ in B^B^ — ► ► i^°7r°)7 [37] is impor- 

tant since it is sensitive to the interference of photons with different polarization, that is, 
photons coming from 7 versus 0' 7 . To illustrate the dependence on the Wilson coeffi- 
cients, we give Sk*j at lowest order (indicated by the superscript (0) for the contributions 
already present in the SM): 

5 *"r = ~TTW ^ " arg( ^° )C ^) ' r = C ' 7,C " ) ■ (fU) 

Here we assume that there is no beyond-the-SM physics in B^ — f^-mixing, and its phase is 
given by the CKM matrix elements. We calculate the exclusive B — ► K*^j decay amplitude 
with QCD factorization following [13] including a s -corrections. The constraints from Sk*^ 
exclude some regions with |r| of order one, unless the CP phases conspire to suppress the 
sine in (6.1), see below. 

The second class of constraints stems from the semileptonic transitions and applies to 

(') ~ 

all Wilson coefficients we consider, C 7 9 10 . The inclusive B — > XJl decays can be predicted 
with high accuracy, in the low-c/ 2 region at the level of < 10% [35], but also the high-g 2 
region is theoretically accessible. As can be seen in Table 3, we utilize the integrated 
branching ratios in the low-g 2 region with q 2 £ [1,6] GeV 2 , B(B — > X s ll)\^^, as well as 
for q 2 > 0.04 GeV 2 , B(B — > X s ll)\ 

[>o.04]- The latter has been experimentally obtained 
by cutting out events with q 2 close to the first and second charmonium resonance, hence 
bears some model-dependence. We use the corresponding theory predictions from [35] and 
[21], respectively. The decay distributions with NP are given in [38]. The treatment of 
uncertainties is as for the B — > X s j branching ratio. 

Concerning the exclusive B — > Kll and B — > K*ll decays, we do not use the branching 
ratios for our model-independent analysis: the constraints are in general weaker than the 
ones from B(B — > X s ll) due to the larger theoretical and experimental uncertainties. A 
particular difficulty with the available exclusive semileptonic decay data is the presence of 
measurements with different dilepton mass cuts, some of which in addition include regions 
where QCDF or SCET does not apply. 

We employ instead early data on the B — > K*ll forward-backward asymmetry from 
Belle and BaBar [4, 5, 6]. While these measurements have large uncertainties, both experi- 
ments strongly support the sign of ^4fb in the high-q 2 region above the second charmonium 
peak to be SM-like. 

A rigorous theory calculation of the exclusive B — ► K*U decays in this kinematical 
region can be facilitated with an operator product expansion in A/Q and m 2 /Q 2 where 
Q = {\^Pi m b} P u t forward in [39]. The leading contribution and also the order m 2 /Q 2 
terms do not introduce new non-perturbative matrix elements beyond naive factorization. 
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Figure 2: Allowed regions of Wilson coefficients in specific NP scenarios after applying the ex- 
perimental constraints from radiative and semileptonic b — > s processes as indicated. The black 
areas are allowed by all constraints. In the left panel we show C' 7 versus Cf F assuming real Wilson 
coefficients. We give the magnitude of C|q P versus its CP phase </>io in the right panel. In both 
plots all other NP Wilson coefficients have been set to zero. For details see text. 



Corrections start to enter at O (a s A./Q). The framework holds at low recoil, (Mb—Mk*) 2 — 
2M#A ^ Q 2 < (Mb — Mk*) 2 , which covers the large dilepton mass region above the ^' 
resonance, q 2 > 14 GeV 2 . 

To leading order in the l/Q-expansion we obtain Afb at low recoil as 



A FB (Q 2 ) oc Re[(Cf (q 2 ) + -^Cf ) C* 10 - (C 9 + -^C 7 ) C 10 ] . (6.2) 

The effective coefficients read as Cf(q 2 ) = C g + (4/3C 1 + C 2 )g{q 2 ) + . . . and Cf = C 7 + . . ., 
where 4/3 C± + C 2 — 0.61 are the dominant SM coefficients. The full expressions including 
the higher order a s -corrections and the QCD penguin contributions are given in [39] and 
are included in our numerical analysis. The lowest order charm loop function is given as 



^ 2 ) = ^ + ^(ln4+"0' ( 6 ' 3 ) 



27 ' 9 V"" q 2 



which agrees with the perturbative quark loop function for massless quarks. Interestingly, 
the dependence on form factors can be factored out in ^4fb (6.2) at this order. We require 
then the sign of (A-fb) integrated over q 2 > 14 GeV 2 to be negative. 

We show the impact of the FCNC constraints on the NP Wilson coefficients for two 
NP scenarios in Figure 2. The areas allowed by all constraints are given in black. We learn 
that the observables (each shown in a different color) yield complementary information, 
and that the SM is allowed, as well as many significantly different NP solutions. 

In the left plot, we entertain NP only in Cj and Cj, and assume further no NP CP 
phases. The regions allowed by B(B — > -XVy), Sk*^ and B(B — > X s ZZ)|r 16 i are shown 
as the green ring, the red cross and the blue half circle, respectively. The impact of 
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Table 4: The ranges of the integrated CP asymmetries (A^ ') for (3min> <Zmax) = (1> 6) GeV 2 after 
applying the experimental constraints at 90% C.L. for the generic NP scenario and those with NP 
in C10, C 10 or C9 only. The upper limits on B(B S — ► /i/z) are also shown. For details see text. 



SK*y is significant. The semileptonic decay excludes in this NP scenario the nipped sign 
solution for the photonic dipole coefficient C7 ~ — C| M ~ 0.31. Note that dimensional 
analysis suggests that power corrections to r of the order C2A / '(3mftCV) ~ 0.1 may induce 
a larger SM contribution to Sk*j than 0(m s /mb) [40]. We show the resulting region in 
the C 7 — C^ p -plane by the dashed lines in the left plot of Figure 2. Because of the present 
experimental situation, however, the inclusion of the power corrections corresponds only to 
a small enlargement of the allowed parameter space. Note also that Ref. [41] estimated the 
non-perturbative SM contributions to be smaller than the ones coming from naive power 
counting. 

In the right plot we allow only for NP in C10, and show the allowed regions in |C]^ P | 
and the CP phase (frio- Fixing the sign of ^4fb (blue) and the semileptonic branching ratios 
B(B — > X S ZZ)|[ 16 ] (green) and B(B — ► -Xy0l[>o.O4] (red) yield orthogonal constraints. An 
upper bound on the magnitude of is obtained with the aid of ^4fb as |C^o P | < 7, 
improving on the bound from the branching ratios alone, \Cy^\ < 10. 

6.2 CP Asymmetries with New Physics 

The dependence of the CP asymmetries on the Wilson coefficients can be seen from the 
analytical (NLO) formulae in Appendix D. We also provide numerical model-independent 
formulae for the B — > K*{-^ Kir)ll branching ratio and CP asymmetries in Appendix F. 
The numerators of ^cp3 4 are sensitive to 67,9 and C 79 whereas the numerators of A® 7 

and Aq probe CV,io and C 710 . The numerators of A^g can be affected by all Wilson 
coefficients considered here. Recall also that y^g are very sensitive to the flipped Wilson 
coefficients since A%$ vanish in the limit C- — > at lowest order in the 1/S-expansion. 
To see directly these features of the T-odd asymmetries, we provide LO formulae: 

A d = 2A D ^(1 - J)Im [(C{? - C 10 )(cf (0) - C 7 )*] , (6.4) 
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Figure 3: The dependence of the integrated T-odd CP asymmetries (A 7 D g g ) for (<Zmin>9max) = 
(1, 6) GeV 2 on NP Wilson coefficients after applying the experimental constraints. In each plot all 
other NP Wilson coefficients have been set to zero. 
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yl 9 = 44 D /3Jlm 



^,(0) w* ^(0)^,/* 4m| eff(o) ^.z* 2mb cff ( ) M w^(o)*^ 

°9 °9 "r" 10 10 i2~W U 7 "i l°7 U 9 — U 7 U 9 / 



^Re(y(°))Im 



2C 7 + ^C£ 



(6.6) 



where for A^,Ag we neglected SM CP violation suppressed by \ u , see Appendix D for 
details. 

We work out the CP asymmetries (A^) with NP by taking into account the experi- 
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mental constraints given in Table 3. We consider scenarios with generic NP, that is, when 
all six NP Wilson coefficients are varied independently, and when varying only one coeffi- 
cient at a time. The asymmetries are integrated over low dilepton masses, q 2 € [1, 6] GeV 2 . 
Theoretical input parameters used are fixed at their central values. 

In Table 4 we show the allowed ranges of the CP asymmetries in various NP scenarios. 
Numerically we find that the CP asymmetries can deviate significantly from their SM 
values, which are doubly Cabibbo-suppressed and below the percent level. As anticipated, 
the T-even CP asymmetries can be enhanced by one order of magnitude up to < 10%. The 
T-odd CP asymmetries A^ g can receive even stronger NP enhancements, up to order one. 

There is also some residual dependence in the on all NP Wilson coefficients from 
the normalization to the CP averaged decay rate. Hence, even though the numerator of 
some CP asymmetries is independent of a particular Wilson coefficient, the asymmetries 
can be modified from their respective SM values given in Table 2. These small effects are 
included in Table 4 whenever they are distinguishable from the SM at la, otherwise called 
SM-like. 

Also the purely leptonic decay B s — > fifi has strong sensitivity to NP contributions 
in 0io and O' 10 since B(B S — ► fifi) oc \C\q — C{ | 2 , see, e.g., [22]. We find a possible 
enhancement of B(B S — ► p,fj,) up to almost an order of magnitude in NP scenarios with 
these coefficients modified, see Table 4. The largest branching ratio, obtained with generic 
NP, is still a factor of two below the current experimental upper bound given in Table 3. 
Furthermore, B(B S — ► Jxji) can be suppressed with respect to the SM by cancellations 
between C\q and C[ . A lower bound exists from data on the decays B — > K^ll, which 
are sensitive to \C\q + C' 10 \ [22]. However, in models containing both and C[ only a 
very weak bound on B(B S — > p,fj,) can be obtained. We conclude that improved data on or 
a discovery of B s — ► jifi decays will have a strong impact on this type of analysis. 

In Figure 3 we show the dependence of the T-odd CP asymmetries (Aj D ^ g ) integrated 

over (q^in, <7max) = GeV 2 on the NP Wilson coefficients as indicated. In the plots all 

other NP Wilson coefficients have been set to zero and the experimental FCNC constraints 

have been applied. The dependence of the asymmetries on the phases is very strong, 

making the CP asymmetries great probes of CP violation beyond the SM. 

In Figure 4 (Aq) is shown against (AP) for different NP scenarios. Both CP asym- 

(') a n 

metries are very sensitive to the phase of C^ , and Aif depends in addition on the dipole 

(') 

coefficients . Correlations of this type can identify the nature of NP. 
7. Summary 

We exploited the full angular analysis in exclusive semileptonic B — > K*(^ Ktt)H and 
B s — > </>(—► K + K~)ll decays as a means of testing the SM and searching for new CP 
phases in b — ► s transitions. From the angular distributions in (2.1) and (2.5) seven CP 
asymmetries A^\ in addition to the one in the decay rate can be accessed. We find that 
the SM predictions, valid for low dilepton masses, have rather large uncertainties ~ 20% 
for A^ 5 , ~ 50% for Acp , Aq , A® s and order one for A$ g , but the tiny magnitude of the CP 
asymmetries < 10~ 2 makes them all ideal to search for a variety of different NP effects. 
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Figure 4: Correlations between the CP asymmetries (Af) and (Aq) in different scenarios of NP: 
generic NP, where all six coefficients are varied, and with NP in Cio and CV only. 

In particular, large NP effects are possible, which survive also the current experimental 
FCNC constraints. We summarize here specific features of the asymmetries: 

- Aj , Ag , Ag are T-odd and can be order one with New Physics. 

- A§ , Aq , Ag , Ag are CP-odd and can be obtained without tagging from dT + dF. 

- ^3,^49 are very sensitive to right-handed currents. 

- A^,Ag, (Aq) can be extracted from a single-differential distribution in (f){0i). 

n (') 

- Af is very sensitive to the phase of the Z-penguins ~ C 1C '. 

Due to the CP-odd feature of some of the asymmetries and the width difference in 
the 5 s -system, further CP asymmetries j\^ D > mix can kg extracted from untagged, time- 
integrated B S ,B S — ► <f>(— > K + K~)ll decays. In the presence of NP with large CP phases 
we expect order one (order ten percent) CP asymmetries related to the T-odd (T-even) 
asymmetries A^ ix (A^ mix ). The difference between the asymmetries originating from 
Bd- and i? s -decays is dominated by the difference in mixing. 

A comparison of the and B s asymmetries hence probes the width difference and 
the mixing phase. One possibility is that the J^ ) ^ m ' ix are order 10 percent, whereas A^ 
are negligible, indicating that there is beyond-the-SM CP violation in AB = 2 processes 
only. Note that a measurement of the ^(^""^ j s probably the only easy way to study CP 
violation in semileptonic rare b — * s transitions with £> s -mesons. 

We conclude that the CP asymmetries from the angular analysis map out precisely 
the CP properties of several Wilson coefficients. The minimal, that is, CKM description 
of CP violation can be disproved. This study can be extended to include also lepton flavor 
dependent effects along the lines of [42]. 
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The prospects for studying rare dimuon modes at the LHC are promising: For an 
integrated luminosity of 2fb~ 1 , i.e., after one nominal year of data taking, a few thousand 
Bd — > K*°fj> + fi~ events are expected at LHCb, allowing a measurement of the branching 
ratio and its CP asymmetry at the percent level [43]. The CP asymmetries proposed here 
require further information on angular distributions, thus higher statistics. A dedicated 
sensitivity study, also taking into account suitable cuts, would be desireable. 



Acknowledgments 



G.H. is happy to thank Frank Kriiger for early discussions on A® and Robert Fleischer 
for careful reading of the manuscript. C.B. would like to thank Sebastian Jager for helpful 
explanations. The work by C.B. and G.H. is supported in part by the Bundesministerium 
fur Bildung und Forschung, Berlin-Bonn. CP. is supported by a grant from the G.I.F., 
the German-Israeli-Foundation for Scientific Research and Development. G.H. is grateful 
to the organizers and participants of the Physics of the Large Hadron Collider-workshop 
and acknowledges the hospitality and great research environment provided by the Kavli 
Institute for Theoretical Physics in Santa Barbara during the final phase of this work. 

A. Angular Coefficients J- a ^ 

Here the functions J- a ^ in the angular distribution (2.1) are given in terms of the transver- 
sity amplitudes Ai_,||,o,t [ 9 ] : 



Jf 



3 j£±® [\A L ± \ 2 + \Af\ 2 + {L^ R)] + ^Re (AlAf + A^Af) }, (A.l) 



= 7W + KT + 



R\2 



Ami 



\A t \ 2 + 2Re{A%A 



L A R* 




h 



3g 
16 

Jjl 
4 

qPi 



3 A 2 



4^2 
3^/2 



\Ai\ 2 + \A\\ 2 + {L^ R) 

\A^\ 2 + (L^R) 
\A L ± \ 2 - \A\\ 2 + {L ^ R) 
Re{A%Af) + (L^R) 
Re{A%Ai*) -{L^R) 



3^2 



Re(AuA 



(L —> R) 



01 



Im(A^Af) -(L^R) 
Im(A%Ai*) + (L^R) 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 
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(A.ll) 



where 



/ 4m? 

A = yi-^, (A.i2) 

and the transversity amplitudes in QCDF can be seen in (C.5). 

B. The effective Hamiltonian 

We use the AB = 1 effective Hamiltonian for b — ► s transitions, e.g., [44, 13] 

n ef f = ~^fV tb V t * s («2 + ^eff) , A u = V ub V: s /V tb V t * s , (B.l) 

where V^- denote CKM matrix elements and we used unitarity to write the basis as 

10 

«2 = c x o\ + c 2 o 2 c + CiOi, n 1 ^ = Ci(o? - or) + c 2 (o 2 c - o?). (b.2) 

i=3 

Here, the 0" 2 denote current-current operators whereas Oj for i = 3, 4, 5, 6 are QCD- 
penguin operators, defined as in [45]. We further take into account the following dipole 
and semileptonic operators 

7 = j^m b [sa^P R b]F^, 0' 7 = j^m b [sa^P L b}F^, 

2 2 
2 2 

^10 = [s^PMh^l], O' 10 = ^—[s^P R b][h^ 5 l], (B.3) 

where Pr/l = (l±75)/2 denote chiral projectors and m b (pL b ) is the MS 6-quark mass at the 
scale fj, b - Since in the SM C[ ~ m s /m b Ci, the chirality flipped operators O 7i910 can only 
compete with 07,9,10 in models beyond the SM. The Wilson coefficients are decomposed 
into their SM and NP parts as Q = Cf M + Cf p and C[ = Of p for i = 7, 9, 10. 

C. Transversity Amplitudes at NLO 

Starting from the K* transversity amplitudes Ai(q 2 ), i = {_L, ||,0, t} in naive factorization 
(see, e.g., [9]), the NLO o^-corrections at large recoil using QCDF [12, 13] can be taken 
into account by the replacements [9, 10] 

(Cf + C'^TM 2 ) - T+, (Cf - C' 7 )TM 2 ) - T~, Cf(q 2 ) - C 9 , (C.l) 
where 

2 F 

Ti ± = ^, T *=Mf^ 7 r = rr+T- (c - 2) 
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The functions ,, can be obtained from the Tj_ 11 given in [12, 13] by substituting Cf with 
C| ff — C\ whereas T± is obtained from 7j_ by replacing Cj S with C| ff + Cy. 

In (C.l), the Tj, i = 1, 2, 3 denote the QCD tensor form factors defined in Appendix E. 
The effective electroweak Hamiltonian employed is given in Appendix B. The effective 
coefficients Cf% and Cg S (q 2 ) have been introduced to absorb 1-loop matrix elements of 
4-quark operators [46]. Here, such contributions to Og are contained in 7^ ± together with 
further corrections beyond naive factorization. We take C| ff and 69,10 at NNLL in the SM 
at the scale fXf,. In the NP scenarios discussed in this work, C' 7 eS equals C' 7 . 

In the framework of QCDF, the functions Tjj, are known at NLO in a s for the SM 
operators and the corresponding chirality flipped operators, see (B.3). The T ± y have the 
following CKM and QCD structure 

T a ± = T a ± V + \ u TM, (C.3) 



r ±(t) = r ±(t),LO , ^± T ±(t). 
a a 4vr a 



NLO 



/Wit) _ n-(u),LO , J2£./t-(u),NLO 
a a An 



where a =_L, ||. At LO in a s (denoted by the superscript (0)) they read 



^±(t),LO _ . 
-'i — sJ 



eff(0) '(0) _^ (o) 

7 7 2m 6 M B 



,7-(u),LO _ f 
^ I - tl 



2m b M B 



y(«)(0) j (c.4) 



r-W' LO = -A, [cf (0) - 4 (0) + ^y(«)l + JT5, ^, (u) ' LO = -6,^yW(°» + HS, 

II L ^JTift II 11 2771,5 

where y((7 2 ) and Y^ u \q 2 ) contain 1-loop contributions of four-quark operators ~ sbqq with 
an imaginary part for q 2 > Am 2 . Since the charm threshold is at the very upper end - if 
not outside - the <7 2 -region where the l/E expansion works and the lighter quarks induce 
either CKM suppressed or penguin contributions, the resulting strong phase is small. In 
(C.4), spectator effects are denoted by HS. At lowest order, these are in -r, (M) ' LO and 



-(t),LO 



. The latter is suppressed by penguin coefficients, whereas the former is non-zero 



only for charged B — > K 11 decays (weak annihilation). At higher order in a s , strong 



phases are further generated in T a 



(i),NLO 



and from spectator interactions [12, 13], which 



have been included in our numerical analysis. The form factors £j_ and £y are discussed in 
Appendix E. 

The transversity amplitudes in the presence of NP Wilson coefficients within QCDF 
and neglecting kinematical terms 2 M\*IM B read as 



A L ^ R = +V2NM B {1 - I) I [(C 9 + C 9 ) T (C 10 + C' w )\ £ ± + 
Ai' R = -V2NM B {\ - s)\ \(C 9 - C 9 ) T (C 10 - C 10 ) 



A 



LtR _ NM 2 (l-s) 



2M K * 



(C 9 - C 9 ) T (C 10 - C( ) - 2m b T- 



2 These formally subleading terms in the l/E expansion are included in the numerical evaluation. 
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NM 2 B {1 - s 



1 (C w — C' w )^ 



(C.5) 



where 



Mi' 



M B 



N 



3-2 10 7r 5 M B ' tb ts 



V, h V*\ 2 sV\f3 l 



1/2 



(C.6) 



and A = M% + M^ t +q 4 - 2(M|Mf » +M B q 2 + M] ( ,q 2 ). Note that A 4 contributes only for 
mi ^ and contains Am, see [12], which represents form factor symmetry breaking QCD 



iL,R 



A^ R for C[ = up to l/E 



corrections. Note that helicity conservation dictates A A 
corrections [47]. 

The dilepton spectrum defined in (2.7) can be written in terms of the transversity 
amplitudes (C.5) as 



dq^ 



\A L ± \ 2 + 



A^\ 2 + \A^\ 2 + (L 



R)+0{m 2 /q 2 ) 



(C.7) 



up to contributions suppressed by the lepton mass. The dependence on the form factors 
stems from £j_ for A^ R , A^ ,R and £y for Aq' R . Since the longitudinal amplitudes Aq ' R are 
enhanced by Mb/Mk*, see (C.5), they imply a stronger dependence of dT/dq 2 on £m than 
on Quantitatively, we find in the SM that for the cuts (<Zmin' 'Zmax) = (1>6) GeV 2 and 
(1,7) GeV 2 the contribution from the longitudinal K* to the total decay rate, 



(\A L ,\ 2 + \A R \ 2 ) 
(dT/dq 2 ) 



(C.8) 



is 0.73±g:?o and 0.7218;??, respectively. 



D. CP Asymmetries and ,4 FB beyond the SM 



Here we give analytical expressions of the CP asymmetries defined in (3.2) and (3.3) in- 
cluding contributions from NP operators (B.3). The asymmetries have been obtained from 
the transversity amplitudes in QCDF, see (C.5), valid in the low dilepton mass region. The 
coefficients C| ff = Cj S,SM + C^ p and Cj are taken into account by T^»- Except for Acp, 
the CP asymmetries are given with their full lepton mass dependence which is confined to 
powers of Neglecting kinematical factors M^„/M B , the CP asymmetries as a function 
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of the dilepton mass can be written as 



A ,8m, r£ M% {l-sf 



Tf\ 2 %7 



+ 



m b \T+\* + \TZ 



+ 



r+ - rr 



-c fl * + 



C* - (<V - -<V) + O (m 2 /g 2 ) , 



(D.l) 



As = .4 



2m b Pi Tt Jm b |T+| 2 - |T_~| 2 7 r - T" 



-Re 



+ 



-Cg* - (5w — ► -<V) 



(D.2) 



= -^l^iRe^ I -^r= s-r- J (C 9 — C 9 )* - 2m b 



2s 



A D = _ A D^ Re 



'T7 X 



14' 



-<V) }, 



c 



^+s5Ljcr -(^^-M 



(D.3) 
(D.4) 



A 6 = A 



^Re(^Icr -^ 



^ = ^Im{(C 10 -C( )^ + ^ -(5 



(D.5) 
(D.6) 



2 ^2 r+(^-) 



£l£|| « 



+ C9C9* + CioCio ~ (<V — ► -<V) j> 



(D.7) 



r+ - rr 



-Cg*- 



T+ + rr 



9 



9°9 ~~ WoWo 



w 



(D.8) 



where (5\v — ► — $w) is short hand notation for conjugating all weak phases. Furthermore, 



r 2 rv 2 

_ Lr F a e 



3 2 • 2 6 ir 5 



\v tb v t * s \ 2 



A 



D 



:\ v tb v ts\ — 



3 2 • 2 6 vr 5 1 



M K *N r 



(D.9) 



where N-p is defined in (3.2). 
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At lowest order in a s , the expressions for the above CP asymmetries simplify by 



q-+,LO _ q ,LO 



^-+,LO , rf- ,LO 



- 2C.' (0) 



7 > 



2cfW + _f_(y(0) + A u y(« 
m 6 



)(o)> 



t-±,lo -r 

-II 
+ s— 



,LO 



= < 



(i _ § )cf (o) + (i + 5)c; (o) 
k (i-5 )( c7 7 efr(o) -c; (o) ) 



(D.10) 
(D.ll) 

(D.12) 



Note that in the SM, or more general, in any model without right-handed contributions to 
the electromagnetic dipole operator, = T±, see Appendix C. 

The lepton forward-backward asymmetry in QCDF is written as 



^4fb 



12AiV 2 Ml(l-5) 2 e 



_L 



dT/dq 2 



(D.13) 



x Re 



m b 



(T+ + T- 



10 



C' + 



c 



10 



E. B — > K* Form Factors at Large Recoil 

The B — ► .ftT* matrix element can be parametrized in terms of seven c/ 2 -dependent QCD 
form factors V, A),i,2 and ^12,3 as 



M B + M x * 
A? 2Mk* 



(E.l) 



(K*( PB - q)\sa^q»(l + l5 )b\B(p B )) = -2i e^e*" p%<f T, 



+ [e*JM 2 B -Ml,)- (e* • q)(2p B - q)^} T 2 + (e* ■ q) 



* ~ Ml- Ml, {2PB ~ q) » 



(E.2) 
T 3 



and 



M B + M K * M B - M K * 

A 3 = 7TT7 A l ^T7 



(E.3) 



2M K * 2M K * 

where e* M denotes the polarization vector of the K* and p B the four momentum of the B 
meson. The QCD form factors obey symmetry relations in the large recoil limit and can be 
expressed at leading order in the 1/E expansion in terms of two universal form factors £j_ 
and £|| [11]. Symmetry breaking corrections at order a s have been calculated using QCDF 
in Ref. [48]. We employ a factorization scheme within QCDF where the are related to 
the V, A 1>2 as [13] 



M B 



M b + M K . 



■v, 



M B + M K » M B - M K * 

S|| = 7^ A l T7 A 2- 



2E 



M B 



(E.4) 
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For the q 2 dependence of the universal form factors we adopt the findings from light 
cone sum rule (LCSR) calculations [25]. Here the q 2 dependence is parametrized as 

n , r 2 



V(q 2 ) 
Ai(q 2 ) 



1 — q 2 /m 2 



+ 



R 



l-q 2 /m 2 it ' 



T2 



] 2 /m 2 



fit 



n 



q 2 /m 2 fit 



+ 



T2 



j 2 /m 2 , y2 



(E.5) 
(E.6) 
(E.7) 



(1 - y /'io fit j 

where the fit parameters r\^,m 2 R and rn 2 it are shown in Table 5. Also given in this table 
are the values of the form factors at q 2 = and the corresponding parametric uncertainties 
within the LCSR approach. We give the uncertainties independent of the Gegenbauer 
moments af'^* and the ones due to af'J^, separately. The relative uncertainty of the form 
factors y(0),^i(0) and A 2 (0) amounts to 8%, 10% and 10% without, and 11%, 12% and 
14% after adding the cl\ : k* ~ see Table 1 for the numerical value - induced uncertainty in 
quadrature, respectively. We use the total relative uncertainty from maximal recoil as an 
estimate for the form factor uncertainties for q 2 > 0. The form factors defined via 
(E.4) are shown as a function of q 2 in Figure 5. Here the bands indicate the uncertainty 
in and £y of 11% and 14%, respectively. 





n 


r 2 m\ [ GeV 2 ] 


m% [ GeV 2 ] 


F(0) 


A F(0) 


A ai F(0) 


V 


0.923 


-0.511 5.32 2 


49.40 


0.411 


0.033 


0.44<5 ai 


A! 




0.290 


40.38 


0.292 


0.028 


0.335 ai 


A 2 


-0.084 


0.342 


52.00 


0.259 


0.027 


0.3l5 ai 



Table 5: The parameters ri^,m 2 R and m^ it describing the q 2 dependence of the form factors V 
and Ai_2 in the LCSR approach [25]. Also shown are the corresponding values of the form factors 
at q 2 = 0, -F(O), their uncertainties independent of the Gegenbauer moment ai.x*, A O -F(0) and the 
uncertainties induced by a\_K» in terms of (5 ai = (ai,x*(l GeV) — 0.1), A ai F(0). 



F. Model-independent CP Asymmetries beyond the SM 



We give numerical formulae for the g 2 -integrated quantities B = t b o (dT/dq 2 ), B = 
t b o (dT/dq 2 ) and Num^^) for Q 2 G I 1 ' 6 ] Gey2 in terms of the NP Wilson coeffi- 
cients Cf p . Here, Num^ip^ denotes the numerators of the CP asymmetries multiplied 
by the 5-meson lifetime such that the normalized CP asymmetries (see (3.5)) are obtained 
from 

Num/A (D) \ 



A 



I 



(F.l) 

B + B y ' 

The dependence of the branching ratios on the NP Wilson coefficients can be written as 



X = Vsm 



1 + NCf T + hReCf p + cTmGf p ) + £ ^Re(Cf p Cf p *)\ for B,B, 

(F.2) 



j>i 
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Figure 5: The universal form factors £j_ and £y in the \ow-q 2 region and their uncertainty indicated 
by the bands. 



whereas the numerators of the T-odd CP asymmetries are parametrized as 
X = X SM [l + (biReCf p + almC? p ) + ^ e i3 -Im(Cf p Cf p *)] for Num (A£ 8 ) . (F.3) 

i j>i 

The numerators of the T-even CP asymmetries read as 



X = Xsm 



1 + {biRjeCf P + CilmCf p ) 1 for Num (A CP>6 ) , Num (Aj 5 ) . (F.4) 



Here, the summations are over i,j = 7, 7', 9, 9', 10, 10' and Xsm denotes the SM prediction 
of the corresponding quantity. Note that for Num (^3,9) we have set Xsm to zero, see 
Section 5, and, hence, the corresponding formulae read as 



X = Y (kReCf p + almCf p ) for Num (A 3 ) , 

i 

X = Y {biReCf p + almCf p ) + ^ Im(Cf p Cf p *) for Num (A £ 

i j>i 



(F.5) 
(F.6) 



The SM predictions Xsm and the coefficients a^, bi, Ci and dij, eij are given in Table 6 and 
Table 7, respectively. We assumed central values for all parameters. 
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Table 6: The SM predictions Xsm and the corresponding coefficients ai, bi and a for i = 
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